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ABSTRACT 

We use the large COSMOS sample of galaxies to study in an internally self-consistent way the 
change in the number densities of quenched early-type galaxies (Q-ETGs) of a given size over the 
redshift interval 0.2 < z < 1 in order to study the claimed size evolution of these galaxies. In a stellar 
mass bin at 10 10 5 < M ga i axy < 10 11 M Q , we see no change in the number density of compact Q-ETGs 
over this redshift range, while in a higher mass bin at > 1O 11 M , where we would expect merging 
to be more significant, we find a small decrease, by ~ 30%. In both mass bins, the increase of the 
median sizes of Q-ETGs with time is primarily caused by the addition to the size function of larger 
and more diffuse Q-ETGs. At all masses, compact Q-ETGs become systematically redder towards 
later epochs, with a (U — V) color difference which is consistent with a passive evolution of their 
stellar populations, indicating that they are a stable population that does not appreciably evolve in 
size. We find furthermore, at all epochs, that the larger Q-ETGs (at least in the lower mass bin) 
have average rest-frame colors that are systematically bluer than those of the more compact Q-ETGs, 
suggesting that the former are indeed younger than the latter. The idea that new, large, Q-ETGs 
are responsible for the observed growth in the median size of the population at a given mass is also 
supported by analysis of the sizes and number of the star-forming galaxies that are expected to be the 
progenitors of the new Q-ETGs over the same period. In the low mass bin, the new Q-ETG appear 
to have ~ 30% smaller half-light radii than their star-forming progenitors. This is likely due to the 
fading of their disks after they cease star- format ion. Comparison with higher redshifts shows that the 
median size of newly- quenched galaxies roughly scales, at constant mass, as (1 + z)^ 1 . We conclude 
that the dominant cause of the size evolution seen in the Q-ETG population is that the average sizes 
and thus stellar densities of individual Q-ETGs roughly scale with the average density of the Universe 
at the time when they were quenched, and that subsequent size changes in individual objects, through 
merging or other processes, are of secondary importance, especially at masses below IO^Mq. 
Subject headings: Galaxies - galaxies: formation - galaxies: evolution. 



1. INTRODUCTION 

The evolution of the median size (i.e., the half-light 
radius ri/ 2 ) of the population of massive (Mcaiaxy > 
1O 1O M ) quenched early-type galaxies (Q-ETGs) at 
given stellar mass has been widely highlighted in recent 
years (e.g., Daddi et al. 2005; Trujillo et al. 2007; Mc- 
Grath et al. 2008; van Dokkum et al. 2008; Cassata et 
al. 2011; Szomoru et al. 2011; Barro et al. 2012; New- 
man et al. 2012; Poggianti et al. 2013, just to cite a 
few). The size of the effect is quite large, with a decrease 
in median r\/2 with increasing redshift oc (1 + z)~ 1- , in 
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coarse terms, this implies that, at a given stellar mass, 
the median half-light radius of Q-ETGs is about a factor 
of ~ 2 — 3 smaller at z <~ 2 than locally, corresponding 
to an increase of over an order of magnitude in the me- 
dian mean stellar density within the half-light radius of 
galaxies. 

The wealth of studies quoted above have used a variety 
imaging data taken from space and from the ground, at 
different wavelengths, and have focused on galaxy popu- 
lations at different redshifts. Quite naturally, there has 
been some debate as to whether obvious observational 
biases might have affected the results, such as the possi- 
ble loss in the noise of outer, low surface brightness parts 
of the galaxies, or the possible effects of color-gradients 
(e.g., Daddi et al. 2005; Mancini et al. 2010). Younger 
stellar populations in the cores of galaxies could result 
in smaller sizes in the rest-frame ultraviolet, where the 
sizes are often measured, than at the longer wavelengths, 
which better sample the stellar mass distribution. Even 
for light profiles measured in the rest-frame optical, the 
strong luminosity evolution of younger stellar popula- 
tions suggests that age gradients within galaxies could 
make the half-light radius of galaxies wavelength depen- 
dent. 

Many of the studies cited above have attempted to 
deal with these uncertainties (e.g., Mancini et al. 2010; 
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Szomoru ct al. 2011), and there is now a reasonable con- 
sensus that there is a real effect to be explained. For 
example, there appear to be no strong color gradients in 
high-z Q-ETG in those cases in which both rest-frame 
UV and optical imaging is available (e.g., Toft et al. 
2007; Guo et al. 2011). Thus, there is now general con- 
sensus that indeed the median size of Q-ETGs is sub- 
stantially smaller at high redshifts, though apparently 
normally sized Q-ETGs coexist with compact ones, es- 
pecially among the most massive galaxies (e.g. Saracco, 
Longhetti & Gargiulo 2010, Mancini et al. 2010; see also 
Onodcra et al. 2010 for a similar conclusion concerning 
the velocity dispersion of Q-ETGs at z = 2). 

With some exceptions (e.g., Valcntinuzzi et al. 2010; 
Cassata et al. 2011; Newman et al. 2012) this trend has 
been often entirely ascribed to the physical growth of 
individual galaxies. Rather than a puff-up mechanism, 
decreasing the central stellar density of Q-ETGs, the fa- 
vored picture has been one in which Q-ETGs maintain 
a nearly constant mass within their innermost few kpes, 
and gradually grow inside-out, building up extended stel- 
lar envclopes/halos around such compact, dense cores 
(e.g., Cimatti et al. 2008; Hopkins et al. 2009, 2010; 
Taylor et al. 2010; Feldmann et al. 2010; Szomoru et 
al. 2011). Accretion of small satellites in minor gas-poor 
mergers has been widely entertained as the leading mech- 
anism to grow these stellar envelops (e.g., Naab et al. 
2009; Hopkins et al. 2009; Feldmann et al. 2010; Nipoti 
et al. 2009; Oser et al. 2012) and thus increase the radius 
of individual high-z compact Q-ETGs, until they reach 
their final z = dimension. 

Yet, this may be only part of the story, and possibly 
a minor one. First, it is now solidly established that 
the population of Q-ETGs has undergone a strong in- 
crease in comoving number density between z ~ 2 and 
the present epoch (e.g., Williams et al. 2009; Ilbert et 
al. 2010; Domingucz Sanchez ct al. 2011). The mass 
functions of different galaxy populations in Ilbert et al. 
(2010), based on high quality photometric redshifts in 
the -2 deg 2 COSMOS field (Scoville et al. 2007), indi- 
cate that the number density of quiescent galaxies has 
increased by a factor of ~ 2 since z ~ 1, and by a factor 
of at least 10 since z ~ 2. These observed number density 
growth factors for Q-ETGs match those expected by ap- 
plying a simple continuity equation to the time-evolution 
of the actively star forming galaxy population (Peng et 
al. 2010,2012). 

In addition, the analysis of SDSS (York et al. 2007) 
mass functions shows that typical passive galaxies with 
M < lO n M0 can only have increased their masses by 
around 20% (and definitely less than 40%) after quench- 
ing (Peng et al. 2010). Analysis of the mass func- 
tions of SDSS central and satellite galaxies in Peng et 
al. (2012) refines these estimates to an average post- 
quenching mass increase of ~ 25% for typical central 
galaxies, and a negligible increase for satellite galaxies. 
These constraints strongly limit the amount of merging 
that may be available to increase the sizes of galaxies. 
This, together with the evolution of the number density, 
can explain why the minor dry merging scenario falls 
somewhat short from quantitatively accounting the ob- 
served size growth since z ~ 2 (e.g., Oser et al. 2012). 

It seems therefore quite likely that the advocated after- 



quenching growth of individual Q-ETGs contributes only 
modestly to the observed secular increase of the me- 
dian size of the Q-ETG population; given the large 
increase of the number density of these systems since 
z ~ 2, it is plausible that another effect may dominate, 
namely quenching of star-formation in actively star form- 
ing galaxies that keeps producing, at later epochs, new 
Q-ETGs with larger size than those of galaxies quenched 
at earlier epochs. The addition, at progressively lower 
redshifts, of progressively larger Q-ETGs will progres- 
sively lower the relative fraction of the more compact 
galaxies relative to the total Q-ETG population, and 
thus produce an upward evolution of the size-mass rela- 
tion. There are good reasons, in an expanding Universe 
that grows structure hierarchically, to entertain the no- 
tion that later-coming Q-ETGs will be larger and thus 
have lower stellar densities than galaxies of similar stellar 
mass that are quenched at earlier epochs (e.g., Coving- 
ton et al. 2011). The apparent disappearance of compact 
Q-ETGs at later epochs may thus be a false reading of 
a reality in which earlier populations of denser Q-ETGs 
remains relatively stable in terms of numbers through 
cosmic time, but become less and less important, in rel- 
ative number, at later and later epochs. 

An important question to answer is thus how the num- 
ber density of compact Q-ETGs evolves from high red- 
shifts all the way down to the local Universe. Searches for 
local analogs to the compact, massive Q-ETGs observed 
at z ~ 2 have been undertaken to answer this question, 
and have given conflicting results. Specifically, for com- 
pact Q-ETGs in massive galaxy clusters, Valentinuzzi et 
al. (2010) have reported evidence for little or no evo- 
lution between z <~ 0.7 and z ~ 0.04 in this population. 
Comparing to the SDSS DR7 (?), other studies have also 
argued for not much evolution in the number density 
of compact Q-ETGs between z ~ 1.5 and the present 
(e.g., Saracco, Longhetti & Gargiulo 2010). Other SDSS 
studies have however reported a drop of at least a fac- 
tor ~ 20 between z ~ 1.6 and z = 0.1 (Cassata et al. 
2011) or even more dramatic than this (e.g., Taylor et 
al. 2010). Also, Szomoru et al. (2011) find that the min- 
imum growth in size required to reconcile the size distri- 
bution of quenched galaxies at z ~ 2 with that of their 
counterparts at z = is a factor ~ 2 smaller than the 
total median size growth observed in the same redshift 
interval. 

Some of the difference between the apparently conflict- 
ing results may in fact stem from compactness having 
been quantitatively defined in different ways by differ- 
ent authors, e.g., either in physical units, or relative to 
the average size-stellar mass relation for local Q-ETGs. 
Quite often, the fraction of compact Q-ETGs is consid- 
ered, as opposed to their number density. Other aspects 
of the analyses need however to come under scrutiny to 
reconcile such widely diverse results. For example, Tay- 
lor et al. (2010), Cassata et al. (2011) and Szomoru et al. 
(2011) use published SDSS sizes as the comparison stan- 
dard at z = 0. This is risky, as the input photometric 
catalogs may have missed compact galaxies through an 
imperfect statistical star-galaxy separation (Scranton et 
al. 2002). Furthermore, as shown by Cibinel et al. (2012) 
on the galaxy sample of the Zurich Environmental Survey 
(ZENS; Carollo et al. 2013a), galaxy sizes smaller than 
the seeing Point Spread Function (PSF) are not reliably 
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recovered from ground-based imaging data. The gener- 
ally poor PSF (FWHM well above ~ 1") of the SDSS im- 
ages casts therefore doubts as to whether the published 
SDSS galaxy catalogues are adequate for this purpose. 
Particularly suggestive is the number density evolution 
of compact galaxies presented by Cassata et al. (2011), 
which, based on the self-consistent analysis of GOODS 
images (Giavalisco et al. 2004) is rather flat from z ~ 2.5 
down to z r~j 0.5, and dramatically drops since z ~ 0.5, 
due to the comparison of the GOODS-based datapoint at 
z = 0.5 with the SDSS point at z = 0. Furthermore, the 
analysis of the number densities in Cassata et al. sums 
up all galaxies above 10 10 M Q , and thus misses possible 
differential evolution with stellar mass. 

The present paper seeks to explore this issue in a care- 
fully controlled fashion by examining, in the redshift 
range 0.2 < z < 1, the evolution at constant size (i.e., 
half-light radius) of the number densities of Q-ETGs, 
and of their plausible star-forming progenitors. We per- 
form our analysis in two bins of stellar mass in which 
the Isi4w < 24 COSMOS sample is complete up to 
z = 1, i.e., 10 10 - 5 - lO n Af and > 1O U M ; these two 
bins straddle across the nearly redshift-invariant char- 
acteristic mass M* ~ 10 n M Q in the Schechter (1976) 
fit to the mass function of galaxies (e.g., Peng et al. 
2011 and references therein), enabling us to search for 
differential effects above and below this mass scale. A 
strength of our analysis is the self-consistent use of data 
from a single survey, i.e., COSMOS, thereby avoiding 
basing our conclusions on comparisons between inhomo- 
geneous samples, and in particular relying on the SDSS 
data for the low redshift reference sample. The COS- 
MOS field is ideal for this study, being unique in having 
both exquisite photometric redshift estimates for a very 
large number of galaxies, based on deep multi-band pho- 
tometry, and high-resolution F814W (I-band) ACS/HST 
images (Koekemoer et al. 2007). While limiting the red- 
shift range between z = 1 and z = 0.2 restricts the evo- 
lutionary lever-arm relative to comparing higher redshift 
samples with SDSS catalogues, our approach has indeed 
the great advantage that the sizes of the galaxies can be 
measured in a uniform way from a single homogeneous 
data set of unparalleled statistical significance. 

In our analysis we use aperture measurements for de- 
termining the sizes of the galaxies from the ACS F814W 
images because of their higher stability relative to model 
fitting approaches when applied to the full morpholog- 
ical diversity of faint high redshift galaxies. We fully 
calibrate however our size measurements against magni- 
tude, size, ellipticity and concentration biases, and show 
that, once both aperture and model-fit measurements are 
so calibrated, they well agree with each other, giving us 
confidence on their robustness. 

In detail, the layout of the paper is as follows. Section[2] 
summarizes the dataset and the basic measurements, ana 
describes the selection criteria for the final galaxy sam- 
ple in detail. Section [3] presents the approach utilized to 
correct sizes and other structural parameters for system- 
atic biases that affect raw measurements as a function 
of galaxy magnitude, size, concentration and ellipticity. 
Section [4] presents the redshift evolution of the number 
densities at constant size- and surface mass density (i.e., 
the size and surface mass density functions) for Q-ETGs, 
and thus our core result, i.e., the constancy of the com- 



pact Q-ETG population and the emergence of a newly- 
quenched population of large ETGs over the z =—> 0.2 
period. Section[5]presents the size- and surface mass den- 
sity functions for star forming galaxies, and compares the 
number densities of the newly- quenched galaxies with 
the number densities of star forming galaxies, of simi- 
lar masses and sizes, that are expected to quench in the 
z = 1 — > 0.2 interval, based on a continuity-equation 
argument (Peng et al. 2010). This Sections also com- 
pares the rest-frame optical colors of compact and large- 
size populations of Q-ETGs, and shows that compact Q- 
ETGs become redder towards later epochs and, at least 
at masses below 10 n Mo, they are also systematically 
redder at any epoch, and thus likely older, than corre- 
sponding large-size Q-ETGs. This reinforces the inter- 
pretation that the latter are the new comers in the Q- 
ETG population, which are responsible for increasing the 
median size of ETGs towards later epochs, without sub- 
stantial increase in size of individual galaxies. In Section 
[6] we conclude. 

Four Appendices present some details of our analy- 
sis. Specifically, Appendices [A] and [B] provide respec- 
tively extra information on the robustness of the mea- 
sured star formation rates and on the reliability of the 
corrections that we apply to the structural/size parame- 
ters; Appendix [C] highlights the general need to correct 
such latter parameters even for estimates based on sur- 
face brightness fitting algorithms which take into account 
the effects of the observational PSF; and Appendix [D] fi- 
nally summarizes the procedure that we follow to derive 
quenching rates using the prescriptions of Peng et al. 
(2010). 

A cosmological model with J7a = 0.7, CIm = 0.3, and 
h = 0.7 is adopted, and magnitudes are quoted in the 
AB system throughout. 

2. THE DATA AND THE BASIC MEASUREMENTS 
2.1. The COSMOS 

We base our stud y on the COSMOS survey dataset 
( |Scoville et"aI1[2007] , so to capitalize on its high quality 
compilation ot multiwavelength imagin g, including Hub- 
ble S pace Telescope (HST) ACS data (Koekemoer et al. 
20071, over a wide-area field. For the present analysis we 
employ the ACS /-band source catalogue of Leauthaud 
et al.| ( |2007[ ) containing 156,748 sources (102,007 of these 
tested to be reliable galaxies) down to a flux limit of 
I814W = 24 mag. The reliability of this catalogue for 
galaxy photometry and morphological analysis was sub- 
sequently improved via extensive visual inspection and 
cleaning to remove artifacts, cosmic rays, and stars, and 
to identify deblending errors, leaving a total of 102,007 
sources flagged as reliable galaxy detections. 

For the purpose of estimating photometric redshifts 
and stellar masses we match the ACS /-band source 
catalo g against the CFHT 
20101) and Subaru 
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j*-band (McCracken et al. 
selected COS MOS Intermediate 
d Broad Band Photometry Catalog ( Ilbert et al.| [2009 



hereafter, the '109 catalog'). This aperture 
matched, photometric database, constructed with an up- 
dated implementation of the source detection procedure 
described in |Capak et al.|2007[ offers a large compilation 
of broad- and narrow- band flux measurements in 3 arcsec 
apertures across 31 bands from UV-optical (u) through 
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to infrared (8.0 /Ltm). As described in Capak et al. (2007) 
the use of PSF-matched, aperture magnitudes allows for 
a single correctional to total flux across all bands for each 
object; estimates of these corrections are pre-compiled in 
the 109 catalog and we adopt these for the present anal- 
ysis. Two further filter-specific corrections to these to- 
tal fluxes are then required prior to photometric redshift 
estimation: (i) a correction against foregro und Galac 



tic dust r eddening, for which we employ the Schlegel et 
al. ( 1998 ) extinction maps with wavelength-de pendent 



Cardelli et al. 



adjustment factor, k\ x E(£? — U), from 
(1989); and (n) a correction against known zero point 
offsets in the COSMOS photometry, for which we em- 
ploy the estimates derived from our photometric redshift 
package ZEBRJ\^\ run in catalog-correction mode (Feld- 
mann et al. 2006)7^ These corrections are consisten t with 
those publish ed fo r the COSMOS pho tometry by |Ilbert| 



et al. 



(2009) and Capak et al. (20071; they are based 



on x -minimization of fitting errors between the best 
fit SED template and observed fluxes of galaxies with 



known redshifts from the zCOSMOS 20k sample (Lilly 
et al. [20071 12009[ ) (see Feldmann et al. 2006) 



For book-keeping purposes we note that our proce- 
dure for matching the HST/ACS /-band source cata- 
log against the 109 multiwavelength photometry cata- 
log (2,017,800 sources) using a 0.6 arcsec tolerance on 
the centroid offset yields a total of 94,908 (93.0%) di- 
rect galaxy matches (i.e., unique galaxy-galaxy associa- 
tions). A further 5267 (5.2%) galaxies were identified 
as sharing their match in 109 with another object in 
the HST/ACS catalog; based on our visual inspection 
of a few hundred such systems drawn randomly from 
the sample, these second matches are typically neigh- 
boring 'junk-sources' or over-deblended fragments of the 
original matched galaxy. We thus treat as successful 
matches the 3614 of these 5267 duplicates for which 
the primary match has a smaller centroid offset against 
its 109 counterpart than the additional (junk) candidate 
match. Conversely, due to the relatively broad PSF of 
the ground-based imaging used in construction of the 
109 catalog, only a small number (16, i.e., ~ 0.02%) of 
galaxies in the HST/ACS-based sample were matched 
against more than one possible 109 counterpart within 
our 0.6 tolerance; once again, after visual inspection of 
these few sources, we adopted as the valid matches the 
sources displaying the smallest offset between the HST 
and 109 centroids. A total of 98,538 galaxies were thus 
deemed successfully matched, leaving only 1816 (1.8%) 
galaxies in the parent HST/ACS catalog unmatched to 
any object in 109. 

One further issue to deal with in the so-obtained galaxy 
sample is the large degree of flux contamination from in- 
terloper objects within the 3" adopted aperture in the 
109 catalog. A total of 13,025 (13%, of our successful 
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9 The four Spitzer mid-IR IRAC bands jSanders et 
are an exception to this rule as it was unfeasible to deg: 
optical data to the much broader Spitzer PSF; rather, the compiled 
IRAC fluxes were measured in fixed apertures of 1.9 arcsec and arc 
corrected to total by dividing out factors of 0.76, 0.74, 0.62, and 
0.58 at 3 .6/xm, 4.5/tm, 5.6/xm, and 8.0/im, respectively; cf. |Ilbert et| 
|al. |2009l 

™ The Zurich Extragalactic Bayesian Redshift An- 
alyzer (ZEBRA) code is available online at our URL: 
www . astro . ethz . ch/ZEBRA . 



matches) were fl agged by Capak et al. (2007) and II- 
bert et al. ( 2009 ) as suffering severe contamination from 
either brighter neighbors or from the diffraction spikes 
of over-exposed stars in at least one of the Subaru Bj, 
Vj, i + , or z + filters. Using such flagged objects may 
introduce errors in our scientific analysis, and hence we 
exclude these systems from our study. The complete- 
ness of our final galaxy sample is consequently ~84% 
(85,513/102,007), contributing a level of uncertainty to 
the absolute normalization of the size- and surface mass 
density £MASS~f unc ti ons & t each epoch comput ed herein 
comparable to that induced by cosmic variance (|Trenti fc| 



Stiavellil2008|) in the COSMOS field ( |Oesch et al.|2010 
We checked that the completeness of our galaxy sample 
does not vary markedly with either size (i.e., half light 
radius) or concentration index, so we do not expect any 
size- or morphology-dependent biases in the presented 
analysis. 

2.2. Photometric Redshifts, Stellar Masses & Star 
Formation Rates 

We estimate photometric redshifts for objects i n our 
matched source c atalog using our ZEBRA code (Feld- 
mann et al.| 2006 ). Calibration of the benchmark SED 



templates (Coleman et al. 1980 Kinney et al. 19961 



employed in this analysis was achieved by comparison 
against a sample of ^20,000 galaxies with secu re spec- 
tro scopic redsh ifts from the zCOSMOS survey (Lilly et 



al. 



2007 2009). Only 236 of the 85,513 input matched 



galaxies were found to be outliers (i.e., a 0.3% failure 
rate). By comparison against the zCOSMOS sample 
at Iguw < 22.5 mag, we estimate a photometric red- 
shift uncertainty of A(z)/(1 + z) ~ 0.007(1 + z) at this 
brightness level. The uncertainty for galaxies down to 
l%\^w = 24 mag was estimated to be 0.012(1 4- z) by 
artificially dimming the photometry of the zCOSMOS 
reference sample. The statistical quality of the ZEBRA 
photo-z's is very similar to that of the COSMOS photo- 
z catalog of Libert et al. (2009). The latter was used 
to further validate the robustness of our results towards 
systematics uncertainties in the photo- z's. 

For each galaxy for which a photometric redshift es- 
timate could be obtained, we further employed a non- 
public extension of ZEBRA (i.e., "ZEBRA+", see Oesch 
et al. 2010) to estimate the corresponding star forma- 
tion rates (SFRs) and stellar masses based on synthetic 
SED fitting to 11 photometric broad-bands, ranging 
from 3832A (u* , CFHT) to 4.5^m (Spitzer IRAC chan- 
nel 2)FM The SED library consists of a comprehensive 
set of star formation history models, i.e., exponentially- 
declining SFRs spanning a range of metallicities from 
0.05 to 2 Zq, decay timescales from r ~ 0.05 to 9 Gyr, 
and ages from 0.01 to 12 Gyr (with a Bayesian prior to 
bound the latter at less than the age of the Universe at 
any given redsh ift). The construction of t his library was 
achieved via the Bruzual & Chariot ( 2003 ) stell ar popula- 
tion synthesis code, adopting a Cha brier IMF ( Chabrier] 
2003 1 . The impact of dust extinction is handled during 
template matching by allowing dust reddening (Calzetti 

11 Note that our stellar masses are defined as the integral of the 
SFR; they are thus about 0.2 dex larger than stellar mass com- 
putations which subtract the mass-return from stellar evolution to 
the interstellar medium. 
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et al. 2000) with the F,(B — V) value treated as a free 
parameter of the fit. Synthetic template matches were 
identified, and a stellar mass successfully derived, for 
all but 1,088 of the 85,277 COSMOS galaxies with pho- 
tometric redshifts (a 1.3% failure rate). Owing to the 
inherent degeneracies in the choice of stellar population 
template and dust extinction model, which dominate the 
error budget in the present analysis (given the thorough 
characterization of the observational SEDs across our 
multiwavclcngth database) we estimate an uncertainty 
of clog a/ ~ 0-20 dex on our model stellar masses. 

For the purposes of separating star-forming from 
quenched galaxies, we adopt a subdivision at an SED- 
fit specific star formation rate (i.e., SFR per unit stellar 
mass; hereafter sSFR) of 10~ n yr -1 ; this corresponds 
closely to the inverse age of the Universe at z ~ 0.3, 
the midpoint of the low redshift bin used in our anal- 
ysis. In Appendix [A] we explain our choice to use star 
formation rates based on SED fits rather than (avail- 
able) IR/UV-derived estimates. In Fig. [T] we show, on 
the rest-frame NUV — J vs. R— J diagram, the distribu- 
tions of quenched and star forming classified galaxies in 
the 0.2 < z < 0.4 and 0.8 < z < 1.0 redshift bins, respec- 
tively; star forming and quenched galaxies are known to 
effectively separate in different re gions in this diagnostic 
plane (e.g. William s et al. 2009, 



Brammer et al. 
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Bundy et al. 2010"[ ). Inspection oi the figure otters con 
fidence that our chosen cut in sSFR well separates star 
forming from quenched galaxies in our sample. 

2.3. Morphological Classification 

The morphological classifications were derived with the 
Zurich Estimator of Structural Types Plus (ZEST+), an 
upgraded version of the ZEST approach described in 
detail in Scarlata et al. (2007a). ZEST+ implements 
well-tested, robust algorithms for measuring a variety of 
non-parametric indices, including concentration, asym- 
metry, dumpiness-, Gini and M20 coefficients for a quan- 
titative structural analysis and morphological classifica- 
tion of faint distant galaxies (see also references in Scar- 
lata et al. 2007a). The new version of the algorithm, 
ZEST+, features several substantial improvements rela- 
tive to ZEST in key computations, including a quality- 
controlled removal of contaminating sources through sub- 
stitution of optimal sky-valued pixels, a more robust 
identification of the galaxy centers, important especially 
for computations of asymmetry and concentration pa- 
rameters, and a more robust calculation of the sources' 
Petrosian radii, unaffected by noise and contaminations 
by nearby sources. ZEST+ also implements a Support 
Vector Machine (SVM) approach to estimate galaxy mor- 
phologies, in addition to the Principal Component Anal- 
ysis (PCA) of the previous ZEST version; the morpho- 
logical classication of the COSMOS sample that we use 
in this work uses the SVM approach. 

Both SVM and PCA method require a training sample 
to guide the morphological classification. The adopted 
training set is classified in three main morphological 
types, i.e., "early- type" (E/S0) galaxies , "disk" (Sa to 
Scd) galaxies, and "very late-type" (Sd/Irr/Pec) galax- 
ies. The galaxies in the training set were carefully se- 
lected as archetypal examples of their classes, well sep- 
arated from non-class members in (at least) concentra- 
tion, asymmetry, and Gini coefficient, and spanning a 
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Fig. 1. — The rest-frame NUV— J vs. R—J color-color distribu- 
tions of quenched (in red) and star forming (in blue) galaxies at 
0.2 < z < 0.4 (top panel) and 0.8 < z < 1.0 (bottom panel) in the 
^814W < 24 COSMOS sample. Galaxies are classified as quenched 
or star forming according to a threshold in sSFR = 10 — 11 yr — 1 . 
This threshold well agrees with the separation in quenched and star 
forming that would be derived using the presented color-color dia- 
gram. To guide the eyes we show as a grey dashed line the saddle 
line between star-forming and quenched samples. Quenched galax- 
ies are largely restricted to [NUV— J — 2.5x (R—J)]> 2.6 mag and 
NUV— J > 4.5 mag. The shown color contours are incremented by 
factors of 2.5, 10, and 25 in number density per bin relative to the 
baseline, which is a factor of 2 lower in the high redshift bin. 



representative range of ellipticities, sizes, and apparent 
mag nitudes (corrected for biases, as discussed in Section 
3.2 1. The morphological classification was performed 
on the reduced HST/ACS Iguw frames. Visual inspec- 
tion of the ZEST+ classified sample, and a quantitative 
inspection based on simulated images, confirms a rela- 
tively small (< 15%) incidence of catastrophic failures in 
the classification, down to and a much higher frequency 
of correctly identified morphologies than in our earlier 
classification attempt using the original version of ZEST 
( |Cameron et al.||2010[ ). 

We emphasize that our main findings are not affected 
by the choice to add a morphological selection to the sam- 
ples of quenched and star forming galaxies. We indeed 
checked that all results stand qualitatively unchanged 
by removing all morphological constraints, with only mi- 
nor quantitative differences which do not affect our main 
conclusions. 



2.4. The Final Galaxy Sample 

The final selection criteria applied to construct the 
master sample of massive galaxies used in the present 
analysis were as follows: (i) ZEBRA maximum likeli- 
hood photometri c re dshift in the interval 0.2 < z < 
1.0 (see Section [2~2|); (n) ZEBRA+ maximum likeli- 
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hoo d ste llar mass greater than 10 M@ (see again Sec 



tion 



2.2 1; (in) ZEST+ morphological type co rres ponding 



to E/SO and Sa-Scd galaxies (see Section 2.3 above) 



(iv) no excessive flux contamination in the groundbased 
CFHT/Subaru imaging from neighboring objects (as in- 
dicated by Capak et al.'s 'bad photometry' flag; see Sec- 
tion [2T|; and (v) no contamination by artifacts, cosmic 
rays, neighboring stars, or from deblending error s in the 
HST/ACS Isiaw imaging (see again Section 2.1). Using 
these selection criteria, and thanks to the excellent com- 
bination of relatively deep ACS images over a ~ 2deg 2 
field, COSMOS returns a final sample of 11,311 (5355 
quenched) galaxies, split in redshifts as follows: 1743 
(921) at 0.2 < z < 0.4, 1751 (833) at 0.4 < z < 0.6, 
3093 (1566) at 0.6 < z < 0.8, and 4724 (2035) at 
0.8 < z < 1.0. The variations of numbers of galaxies 
in the redshift bins highlights that, even COSMOS, is 
not unaffected by cosmic variance. This is a well-known 
fact for this field, and needs attention in order to perform 
studies that involve the evolution of the number densi- 
ties of sources of a given kind. Following the approach 
that we also adopted in Oesch et al. (2010), we correct 
for cosmic variance issues as discussed in Section [4] 

3. GALAXY SIZES: BIASES AND CORRECTION 
FUNCTIONS 

3.1. Raw Size Measurements 

As indicated above, and as also done in the literature 
quoted above, we adopt the half-light radius r 1 / 2 as a 
measure for the size of galaxies in our sample. Whilst 
size measurements based on aperture fluxes are well- 
known to prone to systematic biases in the small-size 
and low surface brightness regim es (e.g. Graham et al. 



2005 Cameron & Driver 



this technique over a proli 



[et~aT1[2007 



20071), we nevertheless favoT 



e- fitting approach (cf. Sargent 
Mancini et al. 2010, Cassata et al. 2011, 
Whitaker et al. 2011) for the present analysis due to 
its high stability — both the stability of its performance 
across the full morphological diversity of the high red- 
shift galaxy population, which becomes increasingly ir- 
regular /clumpy and problematic for mode l fitting codes 
HElmegreen et al.| [2009; Oes ch et al.||2010[ ), and the sta- 
bility (i.e., predictability) of its systematic biases, which 
can th us b e robustly corrected, as we demonstrate in 
Section 3.2 below. Conversely, profile fit-based sizes are 
themselves prone to increasingly unstable behavior (i.e., 
large random errors) in the low surface brightness regime 
( |Haussler et al.|2007[ ) and at faint magnitudes; the results 
of profile-fits are unreliable for typical galaxies fainter 
than Isi4w ~ 23 mag in the COSMOS imaging, whereas 
we want to push our limit one magnitude fainter. 

To measure the (aperture) sizes of the galaxies in 
our sample, we used our cust om-b uilt software package 
package ZEST+ (see Section 2.3). This also measures 
galaxy half light radii within elliptical apertures, pro- 
viding an alternative to the circular aperture measure- 
ments of SExtractor (Bertin & Arnouts 1996), which are 
intrinsically limited in their ability to properly charac- 
terize elongated systems such as inclined disks. Specif- 
ically, given an input measurement of each galaxy's ap- 
parent total flux (we employ the 2.5i?Kron value from 
SExtractor; Kron 1980) and a catalog file highlighting 
the positions of neighboring objects (again from SExtrac- 



tor), ZEST+ estimates the local sky background of the 
galaxy at hand, after replacing the companion galaxies 
segmentation maps with random sky values, and returns 
an estimate for the semi-major axis of the corresponding 
elliptical aperture enclosing half this total flux. 

3.2. Correction Functions for Galaxy Sizes: 
Dependence on Size, Magnitude, Concentration 
and Ellipticity 

It is well established that the brightness of the back- 
ground sky and the blurring effect of the PSF have 
the potential to introduce size and surface brightness- 
dependent biases in both the completeness function of 
galaxy detection and in the re covery of key morpholog- 



ical/structural para meters (cf. Cameron & Driver |2007 
Sargent et al. 20071; hence, the assessment and, il pos- 
sible, correction of any such biases is an essential re- 
quirement for the robust investigation of the evolution 
of galaxy sizes. Despite this general awareness, stud- 
ies of galaxy sizes typically limit themselves at quoting 
global underestimates/overestimates of sizes by a certain 
fractional value in different regimes of brightness and ra- 
dial decline of stellar densities (the latter typically rep- 
resented by Sersic (1963) indices ranging between n = 1 
and n > 4, which more or less bracket the surface bright- 
ness radial declines of exponential disks and elliptical 
galaxies at z = 0). Yet the systematic uncertainties may 
depend also on other parameters, most noticeably the 
ellipticity and size itself of the stellar structure. Here 
we use the approach that we introduced in Cibinel et 
al. (2013) for the ZENS galaxy sample (Carollo et al. 
2013a), and apply the appropriate corrections to mea- 
sured galaxy sizes as a function of galaxy magnitude, 
size itself, concentration, and ellipticity. 

In order to correct for systematic biases in the recovery, 
from COSMOS ACS Iguw imaging, of galaxy flux and 
shape parameters, we performe d > 2,000,000 ar tificial 
galaxy simulations using Galfit (Peng et al. 2002). The 
surface brightness of each model galaxy was defined via 
an elliptical-isophote, Sersic profile with intrinsic (input) 
properties drawn randomly from the following parameter 
space: 14 < Iguw < 25 mag, 0.03 < r 1 / 2 < 6 arcsec, 
< e < 0.95, and 0.2 < n < 8. Each model galaxy 
image was constructed at a pixel scale of 0.03 arcsec, 
matched to our COSMOS/ACS dat a, and convolved with 
a representative ACS Iguw 

PSF flRhodes et alpX^] ). 
The model images were then added with Poisson noise to 
empty regions of sky from the real COSMOS HST/ACS 
images. 

Object extraction was performed on each artificial 
galaxy using an identical approach to that employed in 
constructing the COSMOS Isiaw < 24 mag source cata- 
log. That is, we ran SExtractor first in a "bright source 
mode" (with control parameters: detect_minarea = 
140, detect_thresh = 2.2, and back_size = 400), and 
in the event of null detection a second run was ini- 
tiated in "faint source mode" (with control parame- 
ters: detect_minarea = 18, detect_thresli = 1.0, and 
back_size = 100). The relevant photometric and struc- 
tural parameters output by SExtractor (i.e., the Kron 
magnitude, Kron half light radius, and isophotal ellip- 
ticity) were recorded for each artificial galaxy thus ex- 
tracted. 

The simulation was then advanced to the ZEST+ anal- 
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ysis stage: as noted earlier, we base our analysis on the el- 
liptical aperture galaxy sizes provided by ZEST+, rather 
than circular aperture measurements from SExtractor. 
Hence, as the final stage of our simulation process, each 
artificial galaxy detected above the sky noise was run 
through ZEST+ to quantify its concentration index, el- 
lipticity, and half light radius in a manner consistent with 
our treatment of the real COSMOS HST/ACS sources. 

The suite of artificial galaxy simulations thus compiled 
allows to establish the impact of observational biases 
on the measurement of photometric and structural pa- 
rameters of galaxies in the COSMOS I suw HST/ACS 
catalogue by using the output SExtractor and ZEST+ 
measurements to determine correction functions, for all 
relevant parameters, at each point of a dense grid in 
the four-dimensional space of observed magnitude, size, 
concentration and ellipticity. Specifically, starting from 
regularly-spaced points on a grid in this four-dimensional 
observed space, limited by the above boundaries, we 
selected all models with (output /observed) magnitude 
within 1 mag (Amag — 1), sizes within Ar = 0.015" (or 
Ar = 0.03") for effective radii larger than 0.2"), elliptic- 
ity within Ae = 0.2 and concentration within AC = 0.5 
of the targeted point on the observed grid. We then com- 
puted the median values of the input model parameters 
that had generated these output (observed) values, and 
derived correction vectors connecting the targeted point 
on the output grid with the original point in the input 
grid, i.e., the point with coordinates equal to the me- 
dian of the input values. A correction vector for each 
individual galaxy was then obtained by interpolating the 
correction vectors derived for the grid points. 

The resulting correction functions are illustrated in 
Fig. [2] which shows the grid points in the observed ra- 
dius versus magnitude diagrams, in three bins of con- 
centration (left to right) and ellipticity (top to bottom). 
The length and direction of the arrows give the point- 
by-point median correction vectors. The color of the 
vectors is coded according to the scatter A corr (in ampli- 
tude and/or direction) measured around the median of 
the correction vector (plotted) . The scatter is defined as 
A corr = y/[o-{ri/ 2 )/ < n/ 2 >] 2 + [a(mag) / < mag >] 2 , 
with (7(7-1/2) and < f 1/2 > and a(mag) and < mag > re- 
spectively the 1-a scatter values and medians of the dis- 
tributions of recovered half-light radii and magnitudes at 
the given grid point, respectively. The color of the vec- 
tors varies from green, for A corr < 15%, up to red, for 
A corr > 50%. 

A few well known trends have b een thoroughly com- 
mented on in previ ous studies (e.g. Sargent et al.|[2007j 



|Trujillo et aT 



2007), e.g., the over-estimation of sizes for 
objects with halt light radii < FWHM of the PSF, and 
the under-estimation of sizes (and fluxes) for faint, low 
surface brightness systems with high concentration in- 
dices (Mancini et al. 2010). Inspection of the figure re- 
veals that the required correction functions depend on 
all four parameters (size, magnitude, concentration and 
ellipticity). Applying these corrections to the raw mea- 
surements is important in order not to introduce sys- 
tematic biases in the measurements. In this work we 
make the step of correcting the observed parameters for 
each observed galaxy in our COSMOS/ ACS catalogue, 
by statistically recovering the intrinsic 'true' parameters 



through interpolation between the vectors in the four- 
dimensional calibration grid discussed above. 

In Fig. [3] we show, on the same size-magnitude 
grid points and fixed concentration/ellipticity panels 
("bins"), how systematic measurement errors in concen- 
tration and ellipticity scatter model galaxies out of their 
input/true bins and into incorrect bins of these parame- 
ters. The main point of the figure is to give a "glimpse 
impression" of problematic (yellow to red) regions in 
the four-dimensional parameter space of magnitude, size, 
concentration and ellipticity. In detail, filled green cir- 
cles in the panels indicate grid points in which model 
galaxies simulated with these parameters are more than 
50% of the times correctly recovered within the original 
bin of concentration and ellipticity, and with a negligi- 
ble scatter (AC < 0.3 and Ae < 0.2) relative to the 
true input parameters. Filled red circles indicate grid 
points in which a similarly large fraction of model galax- 
ies simulated with these parameters is recovered within 
the original broad bin of concentration and ellipticity, but 
with a large scatter in concentration index (AC > 0.3) 
and in ellipticity (Ae > 0.2) relative to the true input 
parameters. Squares identify grid points, in a given bin 
of observed ellipticity and concentration, that are dom- 
inated by a fraction of at least 50% (yellow) to > 90% 
(red) of model galaxies which originate in a different bin 
of concentration or ellipticity, and are scattered into the 
observed grid point by measurement errorEpj Finally, 
arrows indicate grid-points from which galaxies that are 
born with those given ellipticity/concentration parame- 
ters 'disappear' from that grid-point, and are observed 
with substantially different ellipticity and/or concentra- 
tion parameters. Specifically, solid arrows indicate points 
of the size-magnitude grid in a given input/true bin of el- 
lipticity and concentration, from which at least 50% of 
model galaxies born in this bin are scattered out from it, 
and end up, due to systematic measurement errors, be- 
ing 'observed' in other (incorrect) bins of ellipticity and 
concentration. Colors from yellow to red indicate an in- 
creasing fraction, above the threshold of 50% required 
to plot an arrow (yellow), of out-scattered model galax- 
ies (with red showing a loss of > 90% of galaxies from 
the grid-point). Empty arrows indicate grid-points from 
which similar fractions of galaxies (> 50%, yellow, up 
to > 90%, red) which are born with those given elliptic- 
ity/concentration values remain within the same elliptic- 
ity/concentration bin, but are recovered with errors in 
concentration AC > 0.3 and ellipticity Ae > 0.2, respec- 
tively. As in Fig. [2] the grey contours show the corrected 
0.2 < z < 1.0, I S uw < 24 mag COSMOS galaxy popula- 
tion(with a change in number density of a factor of two 
between density contours). We discuss in Appendix [B] 
some consistency checks that ensure that the derived cor- 
rection functions return quantities, including sizes, that 
are free from further biases. 

We emphasize that galaxies in some regimes of parame- 
ter space can not be recovered, and thus some systematic 
biases can not be corrected. First, at very low surface 

12 Note that red/green circles can be surrounded by squares, 
indicating that galaxies born with the given grid-point parame- 
ters are well-recovered within the same ellipticity/concentration 
bin, but, to this indigenous population, a population of interloper 
galaxies is added, due to scattering of galaxies into the given grid 
points from other ellipticity/concentration bins. 
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Fig. 2. — Grid of observed half-light radius and magnitude values, obtained with the aperture approach (ZEST+), showing the corrections 
vectors (represented by the arrows) required at each point of the grid to eliminate systematic biases against these same parameters. The 
correction vectors specifically give the extent and direction of the average corrections to sizes and magnitudes required at each specific point 
of the (observed) size-magnitude grid. The nine different panels emphasize the strong influence of profile steepness, and the milder but 
significant influence of profile ellipticity, in shaping the observational biases. In particular, the correction vectors are shown in three bins of 
observed concentration index (left to right: f .25 < c ^ s < 2, 2 < c Q b s < 2.5, and 2.5 < 

c obs < 3.5), and in three bins of observed ellipticity 
top to bottom: < e t, s < 0.3, 0.3 < e ^ a < 0.6 and 0.6 < e b s < 1.0). Green circles indicate points in the observed size-magnitude 
grid in which the corrections are < 0.3 magnitudes and < 25% in radius. The colors of the arrows indicate the degree of uncertainty 
in the amplitude or direction of the arrows, with a scale from green (negligible errors) to red (> 50% error in recovery of the amplitude 
or direction of the correction). Grey open circles indicate regions of the observed parameter space not populated by enough 'observed' 
artificial galaxies to apply a reliable correction (indicating that galaxies which may populate these grid points would be scattered out, by 
observational biasses, into other regions of parameter space; see also Fig. [3] and related discussion in the text). The blue contour lines 
and the grey-shaded contours in each panel show respectively the distributions of galaxies in our final 0.2 < z < 1.0 COSMOS sample 
for uncorrected and corrected measurements (in steps of a factor of two in number density per contour). The dashed line in the top-right 
corner of the panels shows the surface brightness limit for the Iguw < 24 COSMOS sample. 



brightnesses, the original Kron flux under-estimation 
causes objects to fall entirely below our Iguw < 24 
mag catalog selection limit. These regions of classical 
low surface brightness incompleteness are identifiable in 
Figs. [2] and [3j This incompleteness effect has a min- 
imal impact on the observed magnitude-size distribu- 
tions in COSMOS, hampering the observability of only 
the very largest galaxies with magnitudes fainter than 
lauw ~ 23.5 mag, which motivates our initial sample 



selection criteria. Given its nature, this residual bias has 
no impact on the main results that we highlight in this 
paper. Second, at very small sizes, due to PSF broaden- 
ing, galaxies are scattered in regions of low concentration 
and low ellipticity values from bins of high concentration 
and ellipticity; there is thus an intrinsic degeneracy be- 
tween small galaxies which are intrinsically diffuse and 
round, and galaxies which appear so although they are 
in reality more concentrated and elongated than they ap- 
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pear. These residual biases must be kept in mind when 
analyzing galaxy samples; again, they do not affect how- 
ever the results that we present in this paper. 

Finally, in Appendix [C] we show that corrections to 
galactic sizes as a function of magnitude, size itself, con- 
centration and ellipticity are needed not only for aperture 
size measurements, but generally also for sizes that are 
derived from analytical fits to the galaxy surface bright- 
ness distributions. As a showcase we present in this Ap- 
pendix the correction functions, in a similar way as in 
Figures [2] and |3j b ut for size s obtained using the public 
code Gamt ( |Peng et al.|2002 ) . The Galfit analysis in the 
Appendix uses the same set of simulated galaxies which 
were used to correct the ZEST+ size measurements that 
we use in our analysis. We note that the trends and 
strengths of the correction vectors for the Galfit sizes 
differ from those reported in this section for the ZEST+ 
sizes: Galfit does better than ZEST+ in the regime of 
small sizes (i.e., smaller than the PSF), but is less accu- 
rate in recovering sizes of large, low surface brightness 
galaxies. Our main message here is twofold: first, in or- 
der to study galaxy populations spanning a large range of 
sizes and surface brightness, both aperture-based and the 
model-fit size measurements require corrections. Second, 
once these corrections are applied, galaxy sizes derived 
with either of the two approaches are robust and in very 
good agreement with each other, as shown in Appendix 
[C| We therefore stress that the results that we present in 
this paper do not depend on our specific choice for how 
to quantify galaxy sizes. 

4. EVOLUTION IN THE NUMBER DENSITIES OF 
Q-ETGS AT FIXED MASS, SIZE AND SURFACE 
MASS DENSITY 

Fig. [4] shows the number density of Q-ETGs in COS- 
MOS at fixed size (i.e., the size-function &r 1/2 ! lefthand 
panels) and surface mass density Smass (i-e, the Emass - 
function $s MASS ; righthand panels) in four bins of red- 
shift: 0.2 < z < 0.4 (magenta), 0.4 < z < 0.6 (red), 
0.6 < z < 0.8 (orange), and 0.8 < z < 1.0 (grey). 

We restrict our investigation to two fixed intervals of 
stellar mass above the effective completeness limit at 
z = 1 for the spectro-morphological 'early-type galaxy' 
class in the i~si4w < 24 COS MOS catalog (see 



et al. 



Oesch 



2010 Ilbert et al. 20101. In particular, the large 



size of the COSMOS sample enables us to split the to- 
tal sample of massive, > 1O 1O ' 5 M Q-ETGs at these 
redshifts in two distinct mass bins straddling across the 
M GaU 



= 10 U M 







M* value. We thus define a 'low 
mass bin', with 10 10 5 < Afcaiaxy < lO n M (top row in 
Fig. [I]), and a 'high mass bin', with Mcaiaxy > IO^Mq 
(bottom row). We have chosen a linear scale on the y- 
axis to display the size and Smass functions <&; the error 
bars, however, reflect the Bayesian la confidence inter- 
vals (CIs) on the number density in each bin treated as 
a Poisson rate parameter (obtained from inversion of the 
posterior probability distribution for <I>, given the ob- 
served number of objects in that b in and the improper 
uniform prior; see Kraft et al.||1991 for details). 

Note that, in order to correct tor cosmic variance ef- 
fects in the COSMOS field, well known to contain a se- 
vere under-density and a massive cluster at z ~ 0.5 and 
z r-j 0.75, respectively (see, e.g. Scarlata et al. 2007b 



Oesch et al. l 2010 ), we have followed the approach of 
Oesch et al. 1 [2010[ ) and normalized the total number den- 
sity <fr(MGaiaxy > 10 la5 M Q ) of galaxies of all morpho- 
logical types in the 0.4 < z < 0.6 and 0.6 < z < 0.8 red- 
shift bins, and with stellar masses Mcaiaxy > 1O 1O 5 M0, 
against a linear interpolation with redshift of the number 
densities of corresponding galaxies in the 0.2 < z < 0.4 
and 0.8 < z < 1.0 redshift bins. 

To better highlight the global evolution since z = 1 of 
the size- and EMASS-functions, we plot in Fig.[5]the ratios 
between these functions in the lowest (0.2 < z < 0.4) and 
highest (0.8 < z < 1) redshift bins of our analysis. These 
ratios, defined respectively as: 



r — 

1 r 1/2 — 



<fr ri/a (0.2<z<0.4) 
$ ri/2 (0.8<z<l) 



and 



Smass 



<i>s MASS (0.2<z<0.4) 
*s mass (0.8<z<1) ' 



(1) 



(2) 



highlight in a straightforward manner the rate of 
growth in the number densities of Q-ETGs in each bin 
of size and E M ass ■ 

At the highest stellar masses, i.e., in the Mcaiaxy > 
10 11 M e bin of our analysis, we detect a decrease from 
z = 1 to z = 0.2 of about 30%-40% of the smallest and 
densest Q-ETGs, i.e., of high- mass Q-ETGs with half- 
light radii smaller than ~ 2 kpc. Note that in our data 
this is only seen at the 1 — cr- level. Nevertheless this 
result, obtained within a self-consistent dataset, formally 
agrees with the several works which report a decreasing 
number density of compact galaxies with increasing age 
of the Universe (see references above). 

The analysis of our low-mass bin at 10 10 5 < Mcaiaxy < 
10 n MQ, i.e., just below M*, also shows surprises. The 
number density of the 7*1/2 < 2 kpc Q-ETGs remains 
remarkably stable throughout the z = 1 — > 0.2 redshift 
range. At these masses, there is essentially no change in 
the number density of the most compact Q-ETG since 
z ~ 1 down to z ~ 0.2, and thus no evidence for any 
disappearance of compact galaxies due to an increase in 
their individual sizes across this time span. Note that 
we ascertain this constancy self-consistently within the 
COSMOS sample, without recurring to comparisons be- 
tween different studies or datasets. This ensures no spu- 
rious effects introduced by different approaches for com- 
puting the sizes, and by different PSFs, noise properties, 
systematic effects, different depths, and other such po- 
tential complications. 

Similarly striking in Fig. [5] is that, in both mass 
bins, we do observe a strong evolution in the size and 
Emass functions of Q-ETGs, with a marked dependence 
on galaxy size in both mass bins for the population 
growth rates T ri , 2 . This evolution involves however not 
the small-size, compact Q-ETG population, but rather 
the number densities of large-size (and thus low Smass) 
Q-ETGs. In particular, the inspection of Fig. [5] shows 
that, , in both low and high mass bins, a substantial in- 
crease in number density of large-size Q-ETGs has taken 
place since z = 1: the growth rates T ri/2 rise, respec- 
tively in the low mass and high mass bins, from of order 
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Fig. 3. — Grid of true/input half-light radius and magnitude values, obtained with the aperture approach (ZEST+), shown as in Fig. [2] 
split in three bins of observed concentration index (left to right: 1.25 < c ^ s < 2, 2 < c t, s < 2.5, and 2.5 < e |-> s < 3.5), and in three bins oi 
observed ellipticity (top to bottom: < e t, s < 0.3, 0.3 < e t, s < 0.6 and 0.6 < e t, s < 1.0). Filled green circles show grid points in which 
more than 50% of model galaxies born with these given input concentration C and ellipricity e parameters, are correctly recovered with a 
negligible scatter relative to these input parameters (namely, AC < 0.3 and Ae < 0.2). Filled red circles show grid points in which > 50% of 
model galaxies simulated with these parameters is recovered within the original broad bin of concentration and ellipticity (i.e., same panel), 
but with a large scatter in concentration index (AC > 0.3) and in ellipticity (Ae > 0.2) relative to the true input parameters. Squares show 
grid points in which > 50% (yellow, up to > 90%, red) of model galaxies are detected at that grid-point but originate from a different bin 
(panel) of concentration or ellipticity. Solid arrows show grid-points from which at least 50% (yellow, up to > 90%, red) of model galaxies 
born there are scattered out from the given ellipticity /concentration bin (i.e., panel) and are 'observed' in a different panel of ellipticity 
and concentration. Empty arrows show grid-points from which similarly color-codes fractions of galaxies which are born with those given 
cllipticity/concentration values remain within the same ellipticity/conccntration panel bin, but are recovered with errors in concentration 
AC > 0.3 and ellipticity Ae > 0.2, respectively. The direction of the arrows visualizes the global 'scattering direction', due to measurement 
errors, of galaxy models out of that bin. In particular, left /right- ward oriented horizontal arrows indicate recovery-errors in concentration 
only, and specifically scattering towards lower/higher concentrations, respectively. The 45°-inclined left/right- and up/down- ward directed 
arrows indicate recovery errors in both concentration and ellipticity, and precisely towards lower/higher concentrations and lower/higher 
cllipticities, respectively. Remaining symbols in the figure are as in Fig. [2] 
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Fig. 4. — The size (i.e., half-light radius) and surface mass density Smass distributions of quenched early-type, E/S0 galaxies (Q-ETGs) 
in COSMOS, in the redshift interval 0.2 < z < 1. Results are shown for two bins of stellar mass, 10 10 ' 5 < Mgaiaxy < W 11 Mq (top 
row) and MQ a i axy > 1O 11 M0 (bottom row). The redshift interval is split in four bins, i.e., 0.2 < z < 0.4 (magenta), 0.4 < z < 0.6 
(red), 0.6 < z < 0.8 (orange), 0.8 < z < 1.0 (grey). The error bars on each histogram represent only the Baycsian la CIs accounting for 
the Poisson noise in each bin. As noted in Section [2] photometric sample incompleteness and the non-Poissonian component of cosmic 
variance contribute further uncertainties on the order of 15% each to the absolute normalization. Note also that the normalization of the 
0.4 < z < 0.6 and 0.6 < z < 0.8 histograms have been adjusted against cosmic variance in COSMOS, assuming a galaxy number density of 
galaxies with A^Ga'axy > 10 10 ' 5 Mq that linearly increases with redshift, as constrained by the redshift bins 0.2 < z < 0.4 and 0.8 < z < 1.0 
(i.e., excluding the redshift bins 0.4 < z < 0.6 and 0.6 < z < 0.8, which, in the COSMOS field, are known to contain a severe underdensity 
and a massive cluster, respectively). 



unity (i.e., the null growth) and from the negative rate 
of ~ 0.3 at ri/2 < 2 kpc, as reported above, up to ^5-6 
by 7T/2 ^ 5kpc (low mass bin) and 7*1/2 ~ 12 kpc (high 
mass bin). This corresponds to similar growth factors 
from z = ltoz = 0.2 for the population of /ow-Emass 
galaxies, i.e., growth rates Fs MASS ~ 5 from z = 1 to 
z = 0.2 of galaxies with Emass < 9 M kpc -2 . 

The significant increase in number density of large-size 
Q-ETGs with cosmic time between z ~ 1 and z ~ 0.2 im- 
plies that the newly- quenched, large Q-ETGs (hereafter 
NQ-ETGs) lead to a substantial increase in the median 
half-light size for the whole Q-ETG population. The pro- 
portional increase in median size is comparable in the 
high mass and low mass bins, and equal to a factor of 
~1.3 over the redshift span of our analysis. Specifically, 
as shown in Fig. [6] (left panel), the median size of Q- 
ETGs grows from 4.5 ± 0.1 kpc to 6.0±0.3 kpc in the 
high mass bin (solid red line), and from 2.4 ± 0.1 kpc 
to 3.2±0.1 kpc in the low mass bin (dashed red line), 
respectively; these imply a formal growth as approxi- 
mately (1 + z )~ - 70±0 - 15 . Correspondingly, and consis- 
tently with other analyses, there is an increase by a fac- 
tor of ~ 2 — 2.5, between z ~ 1 and z ~ 0.2, in the overall 
number density of M GaJaxy > 10 la5 M Q Q-ETGs. This 
is shown in the right panel of Fig. [6l where we plot the 
number density of Q-ETGs, summed over all sizes, as a 
function of redshift; the result is shown separately for 
our low-mass (10 10 - 5 Af Q < M Galaxy < 10 n M o ; empty 



points) and high-mass (Mcaiaxy > 10 11 M Q ; solid points) 
bins. As discussed in Section [T] these number density 
growth factors for Q-ETGs agree well with the measured 
evolution of the mass function of quiescent galaxies since 
z ~ 1. 

Note that in Figure [6] the increase since z = 1 of the 
number density of Q-ETGs with masses integrated above 
10 n -MQ contrasts the well-known constancy in number 
density of the most massive Q-ETGs, with masses of or- 
der 2-4 x 10 n M o (see also Cimatti et al. 2006; Scarlata 
et al. 2007b; Pozzetti et al. 2010; Ilbert et al. 2013). This 
highlights that the number growth from z = 1 to z = of 
the Q-ETG population at galaxy masses above 10 n M Q 
occurs predominantly within a very narrow mass peak 
around ~ 10 11 M ~ M* . This reinforces the notion 
that the probability for a galaxy to survive quenching 
and keep living as a star forming galaxy drops exponen- 
tially with increasing stellar mass and becomes negligible 
above M* (see Peng et al. 2010). 

Keeping in mind that we observe little or no decrease 
in the number densities of the smallest and most compact 
Q-ETGs since 2 = 1, the key question to answer is what 
is the main cause for the appearance of large Q-ETGs at 
later times - appearance which is largely responsible for 
the measured increase in the median size of the whole 
population of such quenched systems with an early-type 
morphology, whereas the growth of individual galaxies 
has apparently no effect on it. 
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Fig. 5. — The growth factors between 2 = 1 and z = 0.2, r ri/2 and Tj: mass , of the number density of quenched early-type galaxies 
(Q-ETGs) at fixed size (left) and Smass (right)- Results are shown separately for the low mass bin, 10 10 ' 5 < MQ a i axy < 10 11 Mq (top 
row) and for the high mass bin, MQ a i axy > 10 11 Mq (bottom row). The growth factors are defined as the ratios between the size-functions 
(left) and EMASS _mnc tions (right) at 0.2 < z < 0.4, and the corresponding functions at 0.8 < z < 1. The dark-, intermediate-, and 
light-grey vertical bars trace the la, 2a, and 3a Baycsian confidence intervals on these factors, respectively, with the medians marked by 
the magenta+grey datapoints. The coloring of these data points is used as a reminder that they express the ratios between the magenta 
and grey curves in Fig. [" 



"O 

TO 

DC 



X 



"O 



Star-Forming Disks 



10 m5 <M Galaxy <10 11 M s 
M Galaxy >1O 11 M 



Q-ETGs 



~l — 

0.2 



2.5- 



1.5- 



o 
o 

d 1 



0.5- 



Q-ETGs 



O 10 1a5 <M Galaxy <10 11 M e 

• M Ga i axy >10 11 M s 

♦ M Ga , axy >2*10 11 M a 
■ M Ga , axy >4*10 11 M Q 



#5 



0.4 



0.6 
Z 



0.8 



1.0 



0.2 



0.4 



0.6 
Z 



1.0 



Fig. 6. — Left panel: The median size (half light radius) of Q-ETGs (in red) and star-forming disks (in blue), as a function of redshift, 
for the high mass bin (MQ a i axy > 10 11 Mq; solid lines) and low mass bin (10 10,5 Mq < M < IO^-Mq; dashed lines). Right panel: The 
number density of Q-ETGs as a function of redshift in the low- (empty circles) and high-mass (solid circles) bin. Error bars indicate la 
uncertainties on the median values. We also plot the number density evolution of 2 X 1O 11 M0 (diamonds) and 4 X 10 Mq (squares) to 
show that most of the growth in our high mass bin occurs peaked at ~ 1O 11 M0 ~ M* . 



5. DISCUSSION: THE EMERGENCE OF LARGE 
Q-ETGS AT LATE EPOCHS 

We can summarize our findings above with two simple 
statements: (i) the number density of compact < 2 kpc 
Q-ETG remains stable since z = 1, with at most a very 
modest decrease with decreasing redshift of order 30% 
for the most massive > lO n M0 Q-ETGs; and (ii) most 
of the growth of the median size of the > 1O 1O ' 5 M0 Q- 



ETG population is due to the appearance at later times 
of large(r) Q-ETGs. 

On the basis of the current evidence, it is possible of 
course that the Q-ETG population grows in size with 
cosmic time in a coordinated manner, so as to keep con- 
stant the number density of the smallest /densest galax- 
ies, which should be created at the same rate at which 
they would be shifted to the higher size bins. While 
it seems somewhat contrived that such conspiracy may 
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TABLE 1 

Median rest-frame (U — V) colors of compact (low mass bin: < 2 kpc; 

HIGH MASS BIN: < 2.5 KPC) AND LARGE (> 4 KPC) Q-ETGS 



rcdshift 10 10 ' 5 < M/Mq < 10 11 M/Mq > 10 11 

n/ 2 < 2 kpc r 1/2 > 4 kpc r 1/2 < 2.5 kpc r 1/2 >4 kpc 

0.2 < z < 0.4 1.57 ±0.01 1.49 ±0.01 1.56 ±0.04 1.55 ± 0.01 

0.4<2<0.6 1.49 ±0.01 1.45 ±0.01 1.52 ±0.04 1.48 ±0.01 

0.6 < z < 0.8 1.44 ±0.01 1.41 ±0.01 1.53 ±0.02 1.53 ±0.01 

0.8 < z < 1.0 1.42 ±0.01 1.41 ±0.01 1.42 ±0.02 1.50 ±0.01 



TABLE 2 

Median stellar masses within the corresponding mass bins for the 
samples of Table [T] of compact and large Q-ETGs 



10 10 - 5 < M/Mq < 10 11 M/M Q > 10 11 

redshift r i/2 < 2 kpc rj/j > 4 kpc < 2.5 kpc )*i/ 2 > 4 kpc 



0.2 < z < 0.4 
0.4 < 2 < 0.6 
0.6 < 2 < 0.8 
0.8 < z < 1.0 



10.69 ±0.03 
10.71 ± 0.02 

10.69 ±0.01 

10.70 ±0.01 



10.83 ±0.02 

10.82 ± 0.02 

10.83 ±0.03 
10.89 ±0.02 



11.06 ±0.04 
11.11 ± 0.04 

11.07 ± 0.02 
11.07 ± 0.01 



11.34 ± 0.02 
11.32 ± 0.02 
11.36 ± 0.02 
11.32 ± 0.01 



take place, we note that in this case the most compact 
galaxies would be the last that have been quenched, and 
thus they should be systematically younger than aver- 
age. Curr ent analyses of the s tellar popul ation ages in 
the z ~ 2 ([Saracco et al.pOll ) and z = (van der Wei 



et al. 2009f Q-ETG populations find hints for an effect 



in the opposite direction, i.e., the most compact galaxies 
look older, not younger, than their larger counterparts. 

5.1. The rest-frame (U — V) colors of compact and 
large Q-ETGs 

We search for possible trends in stellar population ages 
of compact and large Q-ETGs in our own data: in Fig. [7] 
(top), we plot the distributions of the rest-frame (U — V) 
colors of the low (two leftmost columns) and high (two 
rightmost columns) mass bins, split in two subsamplcs 
with sizes, respectively, < 2 kpc (for the low mass galax- 
ies, and < 2.5 kpc, for the high mass galaxies) and > 4 
kpc. The median colors of both populations in each red- 
shift bin are listed in Table [TJ 

It is immediately evident that there is a trend in the 
same direction as the previous studies at lower and higher 
redshifts. The colors of the most compact Q-ETGs be- 
come on average increasingly redder towards lower red- 
shifts. Specifically, the color difference between compact 
Q-ETGs at z = 1 and at z = 0.2 is A(U - V) ~ 0.16, 
which is in very god agreement with the expected color 
change for a single stellar population that formed at 
z ~ 2 and passively evolves between these two epochs. 
Furthermore, at low masses, compact Q-ETGs are sys- 
tematically redder than their large counterparts at simi- 
lar masses, strengthening the interpretation that the for- 
mer are older than the latter. 

We test that variations in stellar mass within the for- 
mal mass bins do not impact the conclusions above con- 
cerning the relative average ages for the compact and 
large Q-ETG populations. Also in Fig. [7] (bottom) we 
show the distributions of stellar masses, within each red- 



shift bin and formal bin of stellar mass, for the same 
compact and large samples of Q-ETGs of the color anal- 
ysis. The median stellar masses of these distributions 
are listed in Table [2] In each formal bin of mass, the 
median stellar mass of the compact population is very 
constant with redshift. There is understandably a trend 
for the compact Q-ETG sample to have, within each for- 
mal mass bin, a slightly smaller median mass than the 
large Q-ETG population. This goes however in the di- 
rection of decreasing, if anything, the color difference be- 
tween old compact versus young large Q-ETGs. We thus 
conclude that the color difference that we have detected 
between compact and large Q-ETGs in the low mass bin 
is a genuine stellar population effect, indeed consistent 
with a younger average age of large relative to compact 
Q-ETGs. In the high mass bin the two populations, com- 
pact and large, have much more similar colors. At these 
high masses mergers play a much larger role. We will 
further comment on this below. 

In the light of these results, a more plausible interpre- 
tation is that a static rather than dynamic equilibrium 
holds for the number density of < 2 kpc Q-ETGs (i.e., 
the population of compact Q-ETGs remains virtually un- 
changed between z = 1 and z = 0.2, without either cre- 
ation of new compact Q-ETGs or growth of their indi- 
vidual sizes over this time period) . The emergence at low 
redshifts of the population of large and diffuse Q-ETGs, 
seen in both our low- and high-mass bins, therefore re- 
sults from larger star forming galaxies experiencing the 
quenching trauma at later times. 

5.2. Comparison of typical sizes of newly- quenched 
galaxies at different redshifts 

We can use the size functions of Q-ETGs of Fig. [4] to 
derive, in both our mass bins, empirical estimates for the 
size functions of the newly-quenched populations over the 
z = 1 — > 0.2 time period, and their median sizes at the 
median redshift of z ~ 0.6. These estimates are readily 
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Fig. 7. — Top: The rest-frame (U — V) color distributions of Q-ETGs in COSMOS, in the redshift interval 0.2 < z < 1. Results are shown 
for two bins of stellar mass, 10 10 ' 5 < MQ a i axy < IO^Mq (two leftmost columns) and M(3 a i axy > 10 11 Mq (two rightmost column). The 
samples are split in the same four redshift intervals of Fig. [4] Each panel compares the distributions of colors for Q-ETGs larger than 4 
kpc and smaller than 2 kpc (or 2.5 kpc, in the low and high mass bin, respectively). Bottom: The corresponding stellar mass distributions 
for the compact and large samples of Q-ETGs of the top figure. 



defined through the difference between the 0.2 < z < 0.4 
and the 0.8 < z < 1.0 size- functions of Q-ETGs of Fig. [4j 
replotted in Fig. [8] as solid magenta and grey histograms, 
respectively (left and right panels refer respectively to 
the low and high mass bins). The size functions of the 
NQ-ETGs are shown as a black hatched histograms; the 
corresponding median half-light radii of the NQ-ETGs 
are shown as black arrows in Fig. [8] These median sizes 



are equal to 4.3±0.4 kpc and 8.7±0.9 kpc in the low and 
high mass bins, respectively. 

Franx et al. (2008) report a median size of ~ 2 kpc for 
Q-ETGs at z = 2.0 in a mass bin equivalent to our lower 
mass bin. These high-redshift Q-ETGs are all likely to 
have been recently quenched at the time that we ob- 
serve them. Within the large uncertainties involved, the 
change in size of the newly-quenched galaxies is therefore 
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Fig. 8. — Empirical reconstruction for the size functions of NQ-ETGs over the z ~ 1 and z ~ 0.2 time period, defined through the 
difference between the 0.2 < z < 0.4 (magenta histograms) and the 0.8 < z < 1.0 (grey histograms) size-functions of Q-ETGs of Fig.H] Top 
and bottom rows are respectively for 10 10 5 < M < 10 11 Mq and MQ a i axy > 10 11 AfQ galaxies. The black arrows show the median sizes 
of these populations of NQ-ETGs over the z = 1 — > 0.2 period (4.3±0.4 kpc and 8.7±0.9 kpc in the low and high mass bin, respectively). 
In both mass bins, the green histograms show the size functions 3>mq ( r i/2) °f star forming discs predicted to be mass-quenched between 
z = 1 and z = 0.2 according to the continuity equation of Peng ct al. (2010; P10). The green arrows show the median sizes of these 
predicted mass-quenched galaxy populations (5.8±0.1 kpc and 8.7±0.1 kpc in the low and high mass bin, respectively). 



consistent with the idea that NQ-ETGs have, at fixed 
mass, a size that scales roughly as (1 + i.e. with 

the cosmic scaling of the sizes of dark matter halos of 
a given mass. Equivalently, the mean stellar density of 
NQ-ETGs appears to scale roughly as the mean density 
of the Universe. 

This idea may then naturally explain the strong con- 
nection between sSFR and the mean stellar density of 
galaxies in the local Universe. Quenched galaxies with 
low sSFR would naturally have higher stellar densi- 
ties simply because they would, generally, have been 
quenched at significantly earlier epochs. 

5.3. Consistency with predictions from a continuity 
equation 

We finally check whether the global number density 
growth factors that we have observed for the Q-ETGs 
between z = 1 and z = 0.2 are also roughly in agree- 
ment with those expected by applying a number conti- 
nuity equation to the time-evolution of the star forming 
galaxy population, assuming that this is quenched with- 
out a substantial increase of stellar mass in the post- 
quenching phase (Peng et al. 20 10). Anchoring their 
analysis to SPSS JYork et al.|2007l ) and zCOSMOS plly 
et al. 2007 20091 data, these authors show strong evi- 
dence that, while 'environment-quenching' dominates at 
lower masses, the dominant process that piles up galax- 
ies at > 1O 1O ' 5 M0 onto the red-and-dead sequence is a 
quenching mechanism linked only to mass and not to en- 
vironment (i.e., their 'mass quenching'; cf. Figure 13 in 
Peng et al. 2010). 

To explore whether mass-quenching of star-forming 
galaxies is quantitatively consistent, in terms of numbers 
and sizes, with the evolution of the size-function of Q- 
ETGs that we have observed, (z) we compute, in our two 
mass bins, the size (and Smass) functions for the star 
forming population in the COSMOS dataset; and (ii) we 
adopt the predictions of Peng et al. (2010) for sudden 
environment-independent mass-quenching of star form- 
ing galaxies, which is the relevant quenching mode at 



the mass scales of our sample. 

We choose to limit this analysis to the star forming 
disk galaxies in our COSMOS sample. We note that, at 
the stellar masses that we are investigating, disks repre- 
sent > 80% of the star forming population; furthermore, 
we also remark that, above 10 10 5 M Q , the population of 
star forming disk galaxies is largely, > 90%, made of 
systems which have a substantial bulge component. In 
quantitative comparisons with previous studies, it is thus 
important to keep in mind that the star forming popula- 
tion that we consider below is essentially a sample of star 
forming disks with a large if not dominant bulge compo- 
nent. We stress again however that we have tested that 
our results remain unchanged if the morphological selec- 
tions are removed from our study. Apart from relatively 
small differences in the normalizations of the quenched 
and star forming populations, which largely compensate 
each other, the only marginally detected difference is a 
higher number of large-size, irregular star forming galax- 
ies in the high relative to low redshift bin. This is not 
unexpected, and does not alter any of our conclusions; 
on the other hand, the selection of morphologically well- 
behaved samples ensures a higher reliability in the es- 
timates of galactic sizes relative to disturbed, irregular 
morphologies. 

Fig. [9] shows the size- and Emass functions of star 
forming disk galaxies in the same redshift and mass 
bins as in Fig. H| Again we have chosen here a linear 
scale for $ on the vertical axis, and error bars reflect- 
ing the Bayesian la CIs on the number density in each 
bin treated as a Poisson rate parameter. Note that star- 
forming galaxies do not disappear upon quenching, as 
they are continuously replenished from lower mass bins, 
thanks to star formation. Actually, they keep increas- 
ing their number density, as indeed found for their mass 
function both in Ilbert et al. (2010, 2013) and in the Peng 
et al. (2010) model. 

The median half-light sizes of the low mass (dashed 
curve) and high mass (solid curve) star forming popula- 
tion as a function of redshift are furthermore plotted as 
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Fig. 9. — The size (i.e., half-light radius) and surface mass density EjyiASS distributions of star-forming disk galaxies in COSMOS, in the 
redshift interval 0.2 < z < 1. Results are shown for two bins of stellar mass, 10 10 ' 5 < AfQalaxy < 10 11 Mq (top row) and AfG a i axy > IO^Mq 
(bottom row). As in Fig. [4] the redshift interval is split in four bins, i.e., 0.2 < z < 0.4 (dark blue), 0.4 < z < 0.6 (light blue), 0.6 < z < 0.8 
(cyan), 0.8 < z < 1.0 (grey). The error bars on each histogram represent only the Bayesian la CIs accounting for the Poisson noise in 
each bin; photometric sample incompleteness and (the non-Poissonian component of) cosmic variance contribute further uncertainties on 
the order of 15% each to the absolute normalization (Section 12} . As in Fig. [3] the normalization of the two intermediate redshift bins is 
corrected for the known cosmic variance in the field (see text). 



blue lines in the left panel of Fig. [6] Note that the star 
forming disk galaxy population has typically an average 
median size (Smass) a factor of ~ 1.5 — 2 larger (lower) 
than that of Q-ETGs of similar mass and redshift. Also, 
the median size (Emass) of star forming disk galaxies 
shows, within the errors, a remarkable constancy with 
redshift over the z = 1 to z = 0.2 period. These results 
are in qualitative agreement with other work showing 
that the median sizes of star-forming galaxies are larger 
and evolve much slower than those of quenched galax- 
ies (e.g., Franx et al. 2008). Quantitatively, however, 
these and other authors suggest a stronger evolution of 
the average size of star forming galaxies than what we 
find. It is difficult to compare at these later redshifts 
results based on different samples: stellar masses may be 
systematically different, the nominal mass bins may dif- 
fer from ours, and samples typically heterogeneous, i.e., 
they use different galaxy samples to compare different 
redshifts. For example, Franx et al. use, as several other 
authors, the SDSS sample to set the z = reference, 
with few data points at comparably high masses between 
z = 1 and z = 0. We thus refrain from attempting a di- 
rect comparison between our star forming sample and 
others, and stress instead that, using the self-consistent 
COSMOS sample with uniformly derived sizes and well- 
defined morphological properties as described above, we 
actually detect a lack of evolution in the median size of 
10 > 10 10 ' 5 M Q star forming (mostly bulge-dominated) 
disk galaxies over the z — 1 — > 0.2 redshift period. We 
defer to a follow-up analysis to understand the impact 
of cosmological disk fading and other factors on such ap- 
parent lack of evolution in the median size of disk star 



forming galaxies in our sample. 

We then use Eqns. (1) and (27) of Peng et al. (2010) to 
estimate the fractions of star forming galaxies that are 
predicted to undergo mass-quenching over the z = 1 — > 
0.4. Details of the derivation are described in Appendix 
[Dj For the lower mass bin, 10 10 5 < M Galaxy < lO n M , 
this yields quenching fractions /mq of star forming galax- 
ies in the 0.8 < z < 1.0, 0.6 < z < 0.8 and 0.4 < z < 0.6 
redshift bins of 0.22, 0.16 and 0.13, respectively. At 
-^Galaxy > lO n M0, the fractions of newly-quenched 
galaxies in the same redshift bins are instead Jmq = 
0.57, 0.44 and 0.35, respectively. 

We can thus finally compute the size function 
^MQiMcaiaxy,!^^) of star forming galaxies of a given 
mass MQ a i aX y which should be mass-quenched and trans- 
ferred to the Q-ETG sequence between 2 = 1 and z = 0.2 
by summing up the size functions of star forming galaxies 
of Fig. [9] in the three highest redshift bins of our analy- 
sis, eacn weighted by its own mass-quenching fraction as 
given above, i.e.: 

3 

$MQ(fl/ 2 ) = ^2 fMQ,i x $SF,i{ r \/v) ( 3 ) 

with $>sF,i the size function of star forming galax- 
ies in the z-th redshift bin, and i running through our 
three highest redshift bins. The results are presented 
in Fig. [8] as green histograms; their medians are indi- 
cated by the green arrows, and are equal to 5.8±0.1 kpc 
and 8.7±0.1 kpc in the low and high mass bins, respec- 
tively. If mass-quenching of star forming galaxies of sim- 
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ilar mass is responsible for the appearance of the large- 
size Q-ETG populations observed in both our mass bins, 
and if mass-quenching does not change the sizes of the 
progenitors, these predicted distributions for the newly- 
quenched galaxies should be similar to the corresponding 
empirical black-hatched histograms of Fig. IS] 

Remarkably, comparison between black and green 
histograms readily shows that, at masses Mcalaxy > 
lO n M0, the median size for the predicted mass- 
quenched population of star forming disk galaxies over 
the z = 1 — » 0.2 time span (green) is essentially iden- 
tical to that inferred, over the same redshift range, for 
the newly-quenched population of Q-ETGs in COSMOS 
(black). Also the total number densities of 'predicted' 
and 'observed' newly- quenched galaxies in each individ- 
ual size bin between z = 1 and z = 0.2 are in ex- 
cellent agreement, and thus the number densities inte- 
grated over all sizes (respectively 1.1 ±0.1 x 10 -3 Mpc -3 
and 0.9 ± 0.2 x 10~ 3 Mpc~ 3 ). Such good agreement be- 
tween the predicted and observed size functions of newly- 
quenched galaxies is consistent with the idea that mass- 
quenching of star forming galaxies is indeed driving the 
average size evolution of QETGs. 

At lower masses, 10 10 5 < M Ga ia Xy < 1O U M , there 
is still a very good agreement between the number 
densities integrated over all sizes of predicted and ob- 
served NQ-ETGs (respectively 1.3 ± 0.1 x 10~ 3 Mpc -3 
and 1.4 ± 0.1 x 10~ 3 Mpc~ 3 ). At face value, however, 
the median size predicted for mass-quenched star form- 
ing galaxies is larger by about 30% — 40% relative to 
that derived for NQ-ETGs in the z = 1 -> 0.2 period. 
We suggest that the apparent smaller median sizes of the 
newly-quenched remnants relative to their star forming 
progenitors may result from rapid fading of the disk com- 
ponents of these galaxies in the aftermath of their mass- 
quenching episode. This implies that mass-quenching ef- 
ficiently retains the disk components of the star forming 
progenitors, i.e., it has no direct effect of galaxy mor- 
phologies. A similar result is also supported by analysis 
of the morphological mix of the quenched satellite pop- 
ulation in virialized dark matter halos, which is found 
in ZENS (Carollo et al. 2013a) to remain constant with 
halo-centric distance, despite the large increase of the to- 
tal fraction of quenched satellites towards the halo cen- 
ters (Carollo et al. 2013b, in preparation). 

6. SUMMARY AND CONCLUSIONS 

In the previous sections, we have argued that the ap- 
parent size evolution of the Q-ETG population is largely 
due to the appearance of new, larger, galaxies adding 
to a rather stable pre-existing population of more com- 
pact galaxies which has a more or less constant number 
density. We do not claim that this is the only effect in 
operation, but it is certainly in our view the dominant 
one. 

This conclusion is most strongly supported in the lower 
mass bin of Q-ETGs that we have considered in this pa- 
per, 10 10 - 5 < M Ga i axy < lO n M , which samples the 
bulk of the quenched early-type population. At these 
masses, our data show (i) no evolution in the number 
density of the smallest and densest Q-ETGs between 
z = 1 and z — 0.2, and (ii) evidence of a strong build up 
in the population of their larger siblings over the same 



redshift range. This is consistent with the idea that Q- 
ETGs roughly have a constant size, once formed, and 
that newly formed Q-ETG are being added at larger 
sizes. 

The alternative interpretation that galaxies are indi- 
vidually expanding in size while the new ones are added 
to replace them at the small-end of the size distribu- 
tion runs counter to our finding that (i) the compact 
quenched population shows increasingly redder (U — V) 
colors towards lower redshifts, with an average color dif- 
ference between z — 1 and z — 0.2 that is well-consistent 
with the aging, over this time span, of a passively- 
evolving galaxy population, and (ii) that the rest-frame 
ill — V) colors of the larger Q-ETG are systematically 
bluer (as also suggested at the higher and lower redshifts 
respectively by Saracco et al. 2011 and van der Wei et al. 
2009), pointing at a more recent quenching than the cor- 
responding compact population. Together this evidence 
excludes the possibility that the stable number density 
between z — 1 and z = 0.2 of compact Q-ETGs could 
be due to a balance between the formation rate of new 
compact Q-ETGs and their depletion rate due to a size 
growth out of the compact bin, and indicates that the 
newly quenched galaxies are being added at the larger 
end of the size distribution. 

The growth of the median radius of the quenched early- 
type population has usually been interpreted in terms 
of the evolution of individual galaxies (e.g., Oser et al 
2012), without considering the changes in the number 
density that, as argued here, indicate that the effect is 
driven by the addition of newly-quenched larger galax- 
ies at later epochs. The most popular mechanism has 
been size-increase through minor merging. It is clear 
that some such individual size evolution must have hap- 
pened (although our data show that it is not the domi- 
nant driver of the size evolution of the population) . Not 
least, our own analysis has indicated a modest decrease 
in the number density of compact rj/2 < 2 kpc galaxies 
at MGaiaxy > lO n M by - 30% - 40% over the red- 
shift interval z = 1 — > 0.2. These very massive galaxies 
also show weaker color trends; we highlight here, as an 
open issue, to understand the similarity between the av- 
erage rest-frame colors of compact and large quenched 
early-type galaxies with masses above M*. 

A fact to keep in mind is that quenched galaxies that 
are as massive as 10 u Mq and above already by redshift 
z ~ 1 will typically be the central galaxies in haloes that 
have a mass today of > 10 13 Mq. Hydrodynamical simu- 
lations in a A-CDM universe show that, over this period, 
such galaxies grow in stellar mass mostly through minor 
mergers; these mergers produce outer halos and leave the 
central density of galaxies almost constant (e.g., Hopkins 
et al 2009; Feldmann et al. 2010). We also note that 
in the continuity analysis of Peng et al (2010), 10 n M Q 
was identified as a threshold above which post-quenching 
mass increase through mergers was likely to be signifi- 
cant, and below which it was generally unlikely to have 
been. The 30% decrease in density found here, and possi- 
bly the colors, of these very massive quenched early-type 
galaxies are thus consistent with a picture in which a 
fraction of them accrete diffuse halos at later times, as 
also suggested by many other studies (see quoted refer- 
ences). We re-emphasize however that even at this high 
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log(SFR SED ) log(SFR SED ) 

Fig. 10. — Left panel: Fit to the relation between UV-based, dust-corrected SFRs and corresponding values derived from SED fits to 
galaxy photometric points. Right panel: Fit to the relation between the sum of IR-based and dust-uncorrected UV-based SFRs, and SFRs 
derived from SED fits to galaxy photometric points. The good correlation between SFRpjjv,cor and SFRir+uv with SFRsed motivates 
our use of the latter, which is available for the whole Isi4W < 24 COSMOS sample. Note that only one out of 25 points is plotted, for 
clarity of presentation. The contours refer however to the whole sample (with each line referring to an increase in number density of a 
factor of two). 



mass-end the dominant effect in the increase of the aver- 
age radius for the whole population remains the appear- 
ance, at later epochs, of equally massive but larger (in 
half-light size) Q-ETGs. 

The present analysis indicates that the sizes of newly 
quenched ETGs increase with time, and that this is the 
dominant cause of the apparent increase in size of the in- 
tegrated quiescent population. The sizes these of newly 



quenched galaxies appear to scale as roughly (1 + 
suggesting that the mean stellar density within these 
newly quenched galaxies scales as the mean density of the 
Universe at the time of quenching. This idea may unify 
a number of important observational facts, including the 
observed correlation between sSFR and mean stellar den- 
sity in the local Universe. 



APPENDIX 

A. THE CHOICE OF SED-BASED SFRS 

Before adopting as fiducial SFRs those derived from the SED fits, we analyzed in detail SFR estimates derived 
using (a) dust-extinction corrected restframe FUV fluxes, using the Calzetti et al. (2000) extinction curves together 
with the reddening relation E(B - V) = 0.23/3 + 0.49 of Meurer et al. (1999), as well as the Wijesinghe et al. (2011) 
reddening relation E(B — V) = 0.2/3 + 0.33 (with /3 the slope of the UV continuum in a power-law approximation); (b) 
Spitzer/MIPS 24fim detections down to a limiting flux of 80/iJy with Dale & Helou (2002) templates (MIPS catalog- 
courtesy of E. Le Floch); and (c) the combination of IR-based and FUV-based SFRs, the latter not corrected for 
extinction. 

Given the flux limits of the available FUV and Spitzer COSMOS photometry, the FUV- or 24^m-based SFRs result 
in strong and redshift-dependent selection biases towards mass and SFR. The FUV selection picks up 75% of galaxies 
at the redshifts of interest, but it misses the most massive and quiescent systems (median mass of the FUV-detected 
galaxies is log(M/M Q ) = 9.4, compared with a median mass of the galaxies not detected in FUV of log(M/M ) = 9.9). 
Infrared-based SFRs are only available for about 17% of the COSMOS galaxies and pick up only the most massive 
star forming objects (IR median mass < log M/M© = 10.3, compared to 9.3 for galaxies not detected in the MIPS 
images in the 0.2 < z < 1.0 range), with a rapidly increasing star formation rate detection threshold as a function of 
redshift. 

The comparison of the extinction corrected FUV-based SFRs with the SFRs derived by summing the uncorrected- 
FUV plus IR contributions, favors the Meurer et al. reddening relation. Comparison of the FUV SFRs, corrected for 
extinction using this relation, shows a small shift but, overall, an excellent agreement with the SFRs derived from our 
SED fits. A fit of the relation gives: logfS F R F uv,cor) = 1-01 ± 0.01 x \og{SFR SE D) + 0.27 ± 0.01, with a formal 1-cr 
scatter of 0.48 (see left panel of Figure 10 1 . A similar fit to the relation between SED-based SFRs and SFRs obtained by 
summing IR-based and (dust-uncorrecTed) UV-based values gives: log(SFR FUV+IR ) = 1.06 ±0.01 x \og(SFR SED ) + 
0.33 ± 0.01, with a formal 1-cr scatter of 0.45 (see right panel of Figure 10). These excellent agreements validate our 



adoption of the SED-based SFRs for our analysis, which have the virtue to be uniformly derived and available for all 
galaxies in our sample. We stress however that we have tested that our results do not depend on the choice of SFR 
diagnostic adopted for this analysis. 
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B. RELIABILITY OF THE CORRECTION FUNCTIONS FOR GALAXY SIZES: 

CONSISTENCY CHECKS 

Simple analyses and considerations using the uncorrected and corrected (with the corrections of Section[3]) magnitude- 
size distributions shows the importance of accounting for systematic biases in Kron-style size measurements. 

The left panel of figure [TT] shows the distributions of uncorrected and corrected ellipticities of galaxies with concen- 
tration values C < 2, for the entire sample (blue) and the sample limited by r x y 2 < 0.25". Uncorrected measurements 
would lead to an unphysical dependence of sizes of small disk on inclination/viewing angle, as evident from the fig- 
ure. The correction functions restore, from the low ellipticity bin to the high ellipticity bin, a fraction of intrinsically 
small, edge-on disks that were rounded/blurred by the observational PSF, and furthermore recover their intrinsic sizes, 
which were reduced by this PSF blurring. The right panel in the same figure shows the uncorrected and corrected 
size distributions of C > 2.5 galaxies, mostly spheroids and bulge-dominated disks, in two bins of ellipticity, above 
and below e = 0.2. Results are again shown in red and blue for uncorrected and corrected quantities, respectively. A 
significant fraction of intrinsically high ellipticity, C > 2.5 galaxies which our morphological analysis and inspection 
confirm to be primarily edge-on, bulge-dominated disks, are returned by the correction functions to the e cor > 0.2 bin. 
Furthermore, the size distributions of round and inclined systems are in much better agreement after the application 
of the corrections. 



C. 



THE NEED FOR CORRECTING SIZES DERIVED FROM SURFACE BRIGHTNESS FITS: 

A SHOWCASE WITH GALFIT 



As highlighted in Section 3.2 it is not only (Kron, Petrosian, or isophotal) aperture-based measurements which 
require a correction function to aid the recovery of true galaxy magnitudes, sizes, ellipticities, and concentration 
indices from their raw (i.e., observationally biased) counterparts. In this Appendix we illustrate the nature of these 
biases for the case of Sersic (1963) profiler^ modeling with Galfit (Peng et al. 2002) using a process of artificial galaxy 
simulation equivalent to that introduced for the calibration of the SExtractor/ ZEST+ measurements that we adopt 
in out analysis. Specifically, we begin by generating 20,000 artificial galaxy images with input parameters randomly 
drawn within the following limits: 19 < Isi4W < 25 mag, 0.05 < r\/2 < 2.5 arcsec, < e < 0.95, and Sersic index 
0.2 < n < 8. As before each model galaxy was constructed at a pixel scale of 0.03 arcsec and convolved with a 
representative ACS Iguw PSF ( Rhodes et al.|2007 ), and added with Poisson noise to an empty region of sky from the 
real COSMOS HST/ACS imaging. SLxtractor was then run to generate object masks and to estimate starting values 
for use in Galfit. In particular, we adopt the Kron magnitude, the mean isophotal ellipticity and position angle, and 





0.05 0.1 



0.25 0.5 
r-l/2 



Fig. 11. — Left panel: The uncorrected (dashed line) and corrected (solid line) distribution of ellipticities for C < 2 disk galaxies in 
COSMOS. Blue is for the whole C < 2 sample, red is for the subsample of rj/2 < 0.25" galaxies. Right panel: The uncorrected (dashed) 
and corrected (solid) sizes distributions of C > 2.5 disk galaxies in COSMOS. All histograms are normalized to the total number of galaxies 
in the plotted samples. Both panels show that systematics in the raw size measurements are substantially reduced after applying the 
corrections of Section [3] 



13 The well-known Ssersic radial profile is defined as /i(r) = 
fioexp[— (r/o) 1 /"], where ^i{r) is the intensity at radius r and fio is 
the intensity at the galaxy center, i.e., in r = 0. a is the scalelength, 
i.e., the radius at which the intensity drops by a factor e . The 



parameter n is the 'Sersic index', and controls the degree of central 
concentration, i.e., the shape of the profile. 
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the Kron half light radius as inital guesses for their model counterparts — and all fits were begun with an initial guess 
of n — 2.5 for the Sersic index. After running Galfit the output parameters of all convergent fits (i.e., those which 
did not end in a catastrophic error causing the program to cr ash) wer e sa ved for comparison against our input model 
parameters. The results of this analysis are presented in Figs. 



12] and [13] (analogs to Figs. [2] and [3| respectively), 
lations shows two important c 



Examination of the 'correction vectors' diagram for the Galht simulations shows two important differences relative 
to the analog diagram for the SExtractor/ ZEST+ case shown in Fig. [2} (1) not surprisingly, the incorporation of 
the observational PSF into the fitting algorithm Galfit does indeed outperform SExtractor/Zi?5'T+ in size recovery 
for the most compact galaxies; however, (11) in the low surface brightness regime, the Galfit-based measurement are 
substantially affected by systematic offsets from their true (input) values — with a clear trend towards under-estimation 
of both apparent magnitudes and sizes. Note also that, as indicated by the color coding of the error vectors in Figs. 



12 and 13 the random errors in size recovery of high concentration/n-Sersic systems with Galfit are worse in the low 
surface brightness regime than they are for the more stable Kron approach (as implemented in SExtractor/ZEST+). 
Corrections to model-based measurements may be important especially at high redshifts, where faint disk components 
may populate the low surface brightness regime of the correction diagrams. 

Once the Galfit and Kron-based size measurements are properly corrected for, as we do in this paper, t hen both 



estimates obtained with either SExtiax^oi / ZEST+ or Galfit are in excellent agreement, as shown in Fig. 14 The 
figure shows the inherent consistency of our correction functions through a comparison of the apparent magnitudes 
and sizes of thirty randomly drawn high Scrsic-n COSMOS galaxies, recovered after applying our correction maps to 
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Fig. 12. — The correction grid for observed half-light radius and magnitude values for Igi4w < 24 mag COSMOS galaxies, as already 
presented in Fig. reff2, but this time for size and magnitude measurements obtained through Galfit fits to the galaxies' surface brightness 
distributions. The three bins of concentration index of Fig. p] are here replaced by three bins in Sersic index n. Colors and symbols are as 
in Fig. [2] 
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both the SExtractor/ ZEST '+ and the Galfit measurements. Although both measurement techniques require different 
and, in both cases, non-negligible adjustments against systematic recovery biases in this parameter space, the resulting 
distributions are in very good agreement, with only small random offsets between the two sets of estimates. 

This gives us confidence that the results of this study are robust against the choice of size measurement algorithm 
owing to the use of our correction function approach. 

D. THE QUENCHED FRACTIONS FROM PENG ET AL. 2010 

To compare the Peng et al. predictions with our estimates, we convert the mass and time dependent evolution of 
the sSFR in their Eq. (1) : 

S SFR(M G .,„ y ,„ = 2.5 V' (Dl) 

(with f3 = —0.1, see references above), into the mass and time dependent evolution of the reduced specific star- 
formation rate, rsSFR(MG a i axy ,t), as done in Lilly et al. (2013, submitted to ApJ); this is necessary since, in contrast 
with those authors (who use the 'actual' stellar mass trap ped in long-lived stars in their study), we defined the stellar 



mass as equal to the integral of the SFR (see Section 2.2). We have: 
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Fig. 13. — Analysis of scattering of model galaxies in- and out- of true/input ellipticity and concentration bins, as in F ig. reff3, but this 
time for measurements obtained through Galfit fits to the galaxies' surface brightness distributions. Similarly to Fig. |12| the three bins of 
concentration index of Fig. [3] are here replaced by three bins in Sersic index n. Colors and symbols are as in Fig. [3] 
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Fig. 14. — Comparison between uncorrected (red points) and corrected (blue points) parameters obtai ned with ZEST+ (top row) and 
Galfit (middle row). Corrected parameters are obtained by applying the calibrations discussed in Section |3.2| The shown parameters are, 
from left to right: half-light radius r\i^, concentration C, ellipticity e, and Isi4VK magnitude. For clarity, the plot shows results for only 
1500 of the simulated galaxies. The bottom row shows the comparison between measurements obtained with ZEST+ and Galfit before 
(red) and after (blue) the application of the corrections. The applied calibrations substantially improve the scatter between "true" (input) 
and "observed" values, and lead to a very good agreement between corrected measurements obtained with either the ZEST+ aperture 
approach or the Galfit surface-brightness fitting approach. 



rsSFR(M Galaxy , t) = {l-R)x sSFR, (D2) 

where R is the fraction of the baryonic mass that is converted into stars and is promptly returned to the interstellar 
medium (with the remaining fraction 1 — R remaining locked in long-lived stars). The inverse of the rsSFR sets the 
characteristic timescale for the build-up of the long-lived stellar population. As in Lilly et al. 2013, we assume an 
instantaneous mass return and R = 0.4, as derived from stellar population models of Bruzual & Chariot (2003) 



T he fra ction /mq of transient, just being- quen ched objects at any time t can be derived from Eq. (27) of Peng et 



al. (2010), modified to take into account Eq. D2 above, by computing 



, M gsgCt) x r3 SFR(M Galaxy ,t) x %^ m „, 
/mq(<) = $^(t) ' (D3) 

with &sf (t) the number density of star forming galaxies at the same epoch. We can thus predict the fractions /mq,j 
of newly- quenched galaxies in each redshift bin i of our analysis as: 

Imq, = fM Q \lt m '< = f^™" rsSFR(M Galaxy , t)dt, (D4) 

with Zhigh.i and zi ow ^ the upper and lower limit of the i-th redshift bin. 
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